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Neste segundo artigo sobre sistemas Hamiltonianos, apresentamos o metodo da explosao (blow-up) 
para a determinagao da natureza de pontos fixos (pontos de equilibrio) degenerados. Aplicamos o 
metodo a dois modelos hamiltonianos com um e dois graus de liberdade, respectivamente. Primeira- 
mente, analisamos um sistema formado por um pendulo simples submetido a um torque externo 
constante T. Em seguida, consideramos um sistema formado por um pendulo duplo com segmentos 
de comprimentos e massas iguais, tambem submetidos a torques externos constantes e nao nulos. 
A presenga de pontos de equilibrio degenerados nos casos dos pendulos simples e duplo ocorre para 
certos valores dos torques externos. 

In this second article on Hamiltonian systems, we present the blow-up method for the determi- 
nation of the nature of degenerate fixed points (equilibrium points). We apply the method to two 
hamiltonian models with one and two degrees of freedom, respectively. Firstly we study a system 
formed by a simple pendulum submitted to a constant external torque T. Then we consider a 
system formed by a double pendulum of segments with equal lengths and masses, also submitted 
to non-vanishing constant external torques. The presence of degenerate equilibrium points in both 
cases of simple and double pendulums occurs for some values of the external torques. 
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I. INTRODUQAO 

No estudo da dinamica de sistemas Hamiltonianos [H, [2| , a analise do espago de fase tern como 
ponto de partida a busca de pontos de equilibria (tambem denominados pontos fixos). A importancia de 
tais pontos vem do fato de que estes "organizam" a estrutura das orbitas no espago de fase do sistema (ou, 
em outras palavras, reduzem as possibilidades para esta estrutura). Na teoria dos sistemas dinamicos 
um ponto de equilfbrio (ou, de modo equivalente, a estrutura dos espago de fase nas suas vizinhangas, que 
descreve o comportamento dinamico do sistema nesta rcgiao) pode scr classificado mediante o processo 
de linearizagao [H, que consiste na expansao do sistema de equagoes em serie de Taylor ate a primeira 
ordcm, em torno dos pontos de equilfbrio. No cntanto, quando a matriz dos coeficientes constantes do 
sistema linearizado em torno de um dado ponto de equilfbrio (matriz jacobiana) apresenta um ou mais 
autovalores nulos, o processo de linearizagao nao e suficiente para a caracterizagao do ponto de equilfbrio; 
este entao e dito degenerado |5|. 

O matematico Rene Thorn |6| foi o primeiro a obter uma classificagao para pontos de equilfbrio degene- 
rado. Uma classificagao muito simples e tambem apresentada por Bautin Q. Em 1998, Aranda Q apre- 
senta um metodo para analise de pontos de equilfbrio degenerados em sistemas pianos (bi-dimensionais) . 
Um outro metodo muito interessante e o da explosao (blow-up) Q. Estc metodo tern sido aplicado na 
analise da dinamica de diversos sistemas ffsicos com pontos de equilfbrio degenerados. Por exemplo, 
Bogoyavlensky [Ici ] utiliza o metodo da explosao em sistemas astroffsicos e Monerat [Tlj em modelos 
cosmologicos de Friedmann-Robertson- Walker, com constante cosmologica e poeira, para descrever a 
estrutura das curvas no espago de fase do modelo proximas a singularidade. 

Recentemente Monerat et al. [12] apresentaram um estudo analftico sobre o comportamento da 
dinamica dos pendulos simples c duplo, ambos submctidos a torques externos constantes, na vizinhanga 
dos pontos de equilfbrio destcs sistemas. Observa-se a presenga de pontos de equilfbrio degenerados, 
para determinados valores dos torques externos. Para uma descrigao completa desses sistemas, torna-se 
necessario a classificagao destcs pontos de equilfbrio degenerados, objetivo deste trabalho. Chamamos a 
atengao do leitor para a restrigao de nossa discussao a sistemas hamiltonianos; para o tratamento dos 
demais sistemas, sugerimos a referenda @. 

Neste trabalho, faremos uso do metodo da explosao para descrever a estrutura das curvas em uma 
vizinhanga linear dos pontos de equilfbrio degenerados existentes no espago de fase dos pendulos simples 
e duplo. Na segao UH aprescntamos o metodo da explosao, conforme exposto por Guckenheimer e Holmes 
@, para um sistema hamiltoniano de um grau de liberdade, por questao de simplicidade. Na segao 
IIII1 aplicamos o metodo para o pendulo simples submetido a um torque externo constante, descrito por 
uma fungao de Hamilton de um grau de liberdade. Na secao ITVl estendemos o metodo para um sistema 
hamiltoniano de dois graus de liberdade formado por um pendulo duplo, tambem sujeito a agao de torques 
externos constantes. Na segao El aprescntamos nossos comcntarios finais e conclusoes. 



II. O METODO DA EXPLOSAO (BLOW-UP) 

A ideia geral do metodo da explosao e efetuar uma transformagao singular de coordenadas em torno do 
ponto de equilfbrio degenerado, a fim de expandf-lo cm uma csfera n-dimensional contendo um numero 
finito de pontos de equilfbrio do novo sistema. Por questao de simplicidade, apresentaremos nessa segao 
o metodo da explosao para sistemas hamiltonianos com um grau de liberdade, estendendo-o em scguida 
para sistema sistemas com dois graus dc liberdade. Uma descrigao mais abrangente do metodo pode ser 
vista cm 



13] 



Considcre um sistema descrito por uma hamiltoniana TL = 7i(x,p x ) cuja dinamica e governada pclas 
equagoes 

dt ~dp x ~ I[X ' Pxh dt ~ dx~ 9[X ' Pxh (ij 

em que f(x,p x ) e g(x,p x ) sao fungoes das coordenadas {x, p x }. Suponha que exista no espago de fase 
$ do sistema pelo menos um ponto de equilfbrio V de coordenadas x = x*, p x = p*; ou seja, um ponto 
tal que as dcrivadas de x e p x em relagao ao paramctro temporal t sejam nulas. Lembramos que o ponto 
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de equilibrio em questao estara situado sobre uma superficie de energia no espago de fase, definida pela 
equagao E fixa = H{x*, p*). 

Se o ponto dc equilibrio P(x* , p*) for degenerado (ou seja, se a matriz jacobiana do sistema linearizado 
em torno de P(x*, p%) aprescntar pelo menos um de seus autovalores nulos), o comportamento na vizi- 
nhanga linear desse ponto nao podera ser determinado pelo processo de linearizagao do sistema, conforme 
discutido em fl2| . 

A estrutura das curvas em torno do ponto degenerado pode ser analisada pelo metodo da explosao (blow- 
up) 0. Descrevemos os pontos (x, y) uma vizinhanga $ do ponto de equilibrio degenerado P(x*, p*) 
utilizando coordenadas polares (r, 8), tal que r > 0, < 8 < 2ir e 

x — x* = rcosd (2) 
p x - P* x = r scn #- (3) 

Naturalmentc, o caso r = 0, para qualqucr 8, corrcsponde ao ponto V . Apos uma transformagao no 
parametro temporal, as equagoes ([T]) tomam a forma 

Podcmos associar os pontos da vizinhanga $ a pontos de uma superffcie cilindrica T (vide Fig. [1]), de 
coordenadas (cilindricas) (r,0). Scjam T + , Tq c T_ os subconjuntos dc T tais que r > 0, r = c r < 0, 
respectivamente; e seja $ + = $ — {V}. Cada ponto de $ + corresponde a um ponto de r + , e vice-versa. 
O ponto de equilibrio degenerado V corresponde a r , ou seja, a circunferencia r = sobre T. Nao ha 
pontos em $ correspondentes a pontos de T_. 

Determinamos entao os pontos de equilibrio do sistema (|4]), Qi(r = , 8*) e r , tal que i = 1,2, . . . ,n. 
Estes pontos correspondem as diregoes, em segundo as quais as orbitas se aproximam ou se afastam 
do ponto degenerado V . Caso todos os pontos Qi sejam nao-degenerados, a linearizagao do sistema (Q| 
em torno de cada um destes determinara complctamente a estrutura das orbitas, nas suas respectivas 
vizinhangas lineares, e portanto a estrutura das orbitas na vizinhanga de V . Contudo, se ainda houver 
algum Qi degenerado, realiza-se uma nova explosao em torno deste ponto (introduzindo-se novas coor- 
denadas), e assim por diante. Sugerimos ao leitor interessado a Ref.[9( para maiores detalhes e exemplos 
do metodo da explosao. 

III. O PENDULO SIMPLES 



Vamos aplicar o metodo descrito na segao [TT] ao pendulo simples submetido a um torque externo 
constante de intensidade T = mgL, em que m e a massa da particula, L o comprimento do pendulo e g 
a aceleragao da gravidade local. Tal sistema pode ser descrito por uma fungao de Hamilton na forma 



P 

H = — - mgL cos(<p) - mgLip. (5) 



As equagoes de Hamilton que governam a dinamica desse sistema apresentam um ponto de equilibrio no 
ponto P*(ip*,p*), do espago de fase, tal que ip* = |, p* = 0; este ponto possui energia Ef ixa = — T "3 L7r . 
O processo de linearizagao [l2l | das equagoes de Hamilton em torno da vizinhanga linear do ponto P* 
coloca o sistema ((4]) na forma 

d ^-j(x(t)-xA, (6) 



dt 

em que X(t) = (<p{t), Ptp(t)), X* = ((p*, p*) e J e a matriz jacobiana do sistema para o ponto P 




J= — . (7) 



Os autovalores da matriz jacobiana J sao nulos, indicando que o ponto de equilibrio P* e um ponto de 
equilibrio degenerado. 



4 




FIG. 1: O metodo da explosao (blow-up). Em (a), pontos da vizinhanga <£> de um ponto de equilfbrio degenerado 
V sao descritos por coordenadas polares (r, 9) (vide Eq.(|3}); as equagSes de movimento sao escritas nestas coor- 
denadas (vide. Eq.Q). Por sua vez, (b) mostra que (r, 9) podem ser interpretadas como coordenadas cilindricas 
de pontos sobre uma superffcie cilmdrica V. Sejam F+, To e T_ os subconjuntos de V tais que r > 0, r — e 
r < 0, respectivamente; e seja $+ = $ — {V}. Ha uma correspondencia bium'voca entre os pontos de $+ e F+; 
o ponto de equilfbrio degenerado V corresponde a To (a circunferencia r = sobre F); e os pontos de F_ nao 
correspondem a nenhum ponto de 4>. As equacoes de movimento nas coordenadas (r, 9) terao pontos de equilfbrio 
Qi 6 To. Se todos estes forem nao-degenerados, a linearizagao em torno de cada Qi determinara a natureza das 
orbitas em torno de V ■ Se houver ainda algum Qi degenerado, uma nova explosao em torno deste ponto sera 
realizada, e assim por diante. 



Vamos agora aplicar o metodo apresentado na scgao |TTJ Para isso, de acordo com ((3]), faremos uma 
transformagao dc variaveis 

ip = rcos6; p v = rsen0, (8) 
tal que o sistema de equagoes difercnciais formado pelas equagoes de Hamilton assume a forma 

dr — r cos(#) (m 2 <?L 3 sen (rsen(0)) — m 2 gL 3 — rsen(0)) 
T dt = ^L 2 ' (9) 

dO m 2 gL 3 sen(rsen(9))sen(9) — m 2 gL 3 sen(9) — sen(6*) 2 r + r 
T Yt = ^L 2 ' ( ) 



Dc acordo com o metodo, inicialmcntc consideramos o limite r — > nas equagoes Q c (|10[) ; em seguida, 
dctcrminamos os valores de 9 = 9* que fornecerao as coordenadas do ponto de equilfbrio degenerado nas 
novas coordenadas. O resultado neste caso e 9* = 0. Tal ponto esta associado a uma superffcie de encrgia 
E/ixa = —mgL. Entao, de acordo com o metodo exposto na segao|TIl linearizamos o sistema em torno 
do ponto (r = 0, 9* = 0); a matriz jacobiana nas novas variaveis e 

/ mgL \ 
J= . (11) 

V —mgL 1/mL 2 I 

A natureza de tal ponto dc equilfbrio e determinada pelos autovalores da matriz J: 
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_ 1 + y/1 - Am 4 L e g 2 _ 1 - yj - 4?n 4 L 6 g 
1 2mL 2 ' 2 2mL 2 ' [ ' 



5 



Podemos observar que se to 4 = 4L i g 2 ambos os autovalores sao reais e Ai = A2 > 0, indicando que 
tal ponto apresenta um equilibrio instavel, denominado "no flexionado" (inflected node) [J. A condigao 
to 4 > 4j J} g5 indica que o ponto de equilibrio e do tipo Foco, ou seja, autovalores complexos da forma 
A = a ± ib, em que {a, b} sao mimeros reais e positivos 3. Esse resultado indica que o volume do espago 
de fase nao e conservado. Em sistcmas Hamiltonianos conservatives, com dois graus de liberdade, pode 



ocorrer o aparecimcnto de pontos de equilibrio do tipo Foco-Generalizado conforme mostrado em 14 1. 
Nesses casos cstes pontos representam situagoes de equilibrio instavel. Para m 4 < A] \ g i os autovalores 
sao reais e distintos, configurando um no instavel. 

A solugao geral do sistema linearizado formado pelas equagoes © c (|TU1) em torno do ponto de equilibrio 
e uma superposigao (combinagao linear) das solugoes linearmente independentes: 



m— 1 

em que os A m ' sao os autovetores associados aos autovalores Am da matriz Ji, e os coeficientes Cm sao 
constantes de integragao que dependem das condigoes iniciais escolhidas. 

Em termos das variaveis r, 9 as solugoes [l^ validas numa vizinhanga linear do ponto de equilibrio sao 
da forma 

r{t) = ciA ie Al * + c 2 e X2t (14) 
0(t) = c 3 A 2 e Alt + c 4 e A2i . (15) 

em que Ai e A2 sao dados por (fT2j) . Como ja mencionado, as solugoes (fl"4f c (fl~5|) dependem da relagao 
entre os valores da massa e do comprimento do pendulo. A estrutura das curvas na vizinhanga linear 
dos pontos de equilfbrios para os casos de no instavel (quando m 4 < 4L \ a ), no flexionado (quando 

m 4 = jjjg-r) e foco (quando m 4 > jjj-^) podem ser vistos na Ref.0]. 

IV. O PENDULO DUPLO 

Um pendulo duplo de massas identicas m e segmentos de comprimentos identicos L, sob a agao de 
torques externos constantes de intensidades fti = 2mqL e P2 = rngL, respectivamente, apresenta um 
ponto de equilibrio degenerado no seu espago de fase[12| . Este sistema e descrito por uma fungao de 
Hamilton de dois graus de liberdade na forma 

n, T+^i _C0S (Vl _ ^2)piP2 (-) t { \ r I \ , o™ r r na\ 

H = — — r^r 2771.9^008(^1) - mgLcos((p 2 ) + 2mgLip 1 + mgLtp 2 . (16) 

As equagoes de Hamilton governam a dinamica do sistema; elas constituem um conjunto de quatro 
equagoes diferenciais nao lineares, e possuem um ponto de equilibrio degenerado Pq de coordenadas 

P : (pi = 0, Vl = -it/2, p 2 = 0, <p 2 = -vr/2), (17) 

com uma energia associada E = —^irmgL. Em se tratando de um ponto de equilibrio degenerado, o 
processo de linearizagao do sistema de equagoes em torno desse ponto nao e suficiente para caracteriza- 
lo (i.e., caracterizar o comportamcnto do sistema em sua vizinhanga). Para dcterminarmos a natureza 
deste ponto de equilibrio faremos uso do metodo da explosao, apresentado na segao HIl O pendulo 
duplo possui dois graus de liberdade (1^1,(^2), aos quais associamos os momenta canonicamente conju- 
gados [pi,p 2 ), respectivamente. Portanto, utilizarcmos coordenadas hipercsfericas quadridimensionais, 
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(R, 9, (p, 77), definidas pelas transformagoes 1 de variaveis 

ipi = R sen(0) cos(0) sen(?7) ; pi = R sen(9) sen(cj)) seu(rj), 

<P2 = R cos(#) sen(?7) ; p 2 = R cos(^), 



(18) 



cm que 77 £ [O,tt],0 £ [0, 2ir],(j) £ [0,2it] c R £ [0, 00). Substituindo as cxprcssoes (fl8|) nas equagoes de 
Hamilton do sistema obtido a partir da hamiltoniana (|16j) . obtemos urn novo conjunto de quatro equagoes 
diferenciais nas novas variaveis (R, 9, </>, 77). As equagoes de Hamilton nas variaveis originais pode ser visto 
na Ref. fl2l |; nas novas variaveis, as equagoes tomam a forma 



R w = AOM^); Rf = f 2 (R,e,4>,v) 



R 



dt 



MR,0,<f>,ri)) R$L = U(R,9,<I>,ti). 



(19) 



As coordenadas do ponto de equilibrio degenerado nas novas variaveis sao obtidas anulando o lado 
esquerdo das equagoes (fl"9)) . e em seguida, considcrando o limite R — > das fungoes i £ {1,2,3,4}. 
Temos dessa forma um sistema algebrico de tres equagoes (ja que a primeira equagao e identicamcnte 
nula) , 



—2mgL cos(cj)) [sen(0)sen(r7)] _1 = 0; 
—2mgL sen(^) cos(#) [sen^)] -1 = 0; 
~2mgL cos(rj) sen(#) scn(</>) + mLg sen(?;) = 0, 
cuja solugao fornece quatro pontos de equilibrio, 



(20) 



' Qi : (R = 


0, 


= tt/2, v = 


arctan(2), 9 — 


t/2); 


Q2 ■■ (R = 


0, <j> 


= tt/2, V = 


— arctan(2), 9 


= -tt/2); 


Q 3 : {R = 


0, <f> 


= -7r/2,r) 


= — arctan(2), 


9 = tt/2) ; 


k Qi : (R = 


o,<j> 


= -ir/2, V 


= arctan(2), 9 - 


= -7T/2). 



(21) 



Todos os pontos de equilibrio descritos em (|2"Tj) estao associados a mesma superficie de energia 
e,: = —'SmgL. Ao linearizarmos o sistema de equagoes (fl"9|) em torno de qualquer um dos pontos de 
equilibrio, obtemos matrizes jacobianas identicas, da forma 



Ji 



( -y/EmgL 

VEmgL 

^ImF 1 VSmgL 

V y/lmgL j 



(22) 



em que i £ {1,2,3,4}. O conjunto de autovalores das matrizes jacobianas fornecerao a natureza do 
i-esimo ponto. Obtemos, 



A1.2 = ±VEmgL] A3, 4 = VbmgL. 



(23) 



Uma hiperesfera quadridimcnsional de raio R e o lugar geometrico dos pontos do espago euclidiano 5ft 4 , da forma 
(<pi, pi, 1^2, P2), que guardam a mesma distancia _R da origem, satisfazendo portanto a equagao cartesiana ip 2 + p 2 + 
(^2 + P2 = Esta equagao e a generalizagao do caso bidimensional (circunferencia, x 2 + y 2 = R 2 ) e tridimensional 
(superficie esferica, x 2 + y 2 + z 2 = R 2 ). 
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Este resultado indica |4| que os quatro pontos de equilfbrio correspondem a combinagao de uma sela 
(^1,2) e de um no flexionado (A34); representam, portanto, situagoes de equilfbrio instavel. 

V. CONCLUSAO E COMENTARIOS FINAIS 

Ncstc trabalho, aplicamos o metodo da explosao (blow-up) para a determinacao da natureza de pontos 
dc equilfbrio degenerados no espago de fase de dois modelos hamiltonianos: o pendulo simples e o pendulo 
duplo, ambos sob a acao de torques externos constantes. Conclufmos que todos os pontos degenerados 
existentes nestes dois sistemas descrevem situagoes de equilfbrio instavel. No caso do pendulo simples, 
mostramos que ha apenas um ponto de equilfbrio dcgencrado; este pode scr um no instavel simples ou 
flexionado. dependendo da relacao entre a massa e o comprimento do pendulo. No caso do pendulo duplo, 
apos a explosao observamos a existencia de quatro pontos de equilfbrio instaveis de mesma natureza; cada 
um deles consistem no produto direto de uma sela hipcrbolica (associada a um par de autovalores rcais 
e simctricos) por um no flexionado. Assim, o metodo em questao permite a descrigao das solugoes na 
vizinhanga linear dc cada um dos pontos de equilfbrio degenerados, mostrando ser uma tecnica eficaz na 
analise da estabilidade/instabilidade de tais sistemas. A relativa simplicidade do metodo permite que 
este possa integrar o programa de uma disciplina de Sistemas Dinamicos a ser cursad a ap os a disciplina 
de Mecanica Analftica[15| nos cursos de graduacao em Ffsica, conforme discutido em 
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